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Ferrofluids are strongly paramagnetic liquids. We study the behavior of ferrofiuid droplets 
confined between two parallel plates with a weak applied field parallel to the plates. The droplets 
elongate under the applied field to reduce their demagnetizing energy and reach an equilibrium shape 
^ ^ . where the magnetic forces balance against the surface tension. This elongation varies logarithmically 

^ ^ with aspect ratio of droplet thickness to its original radius, in contrast to the behavior of unconfined 

droplets. Experimental studies of a ferrofiuid/water/surfactant emulsion confirm this prediction. 



< 



^ I. INTRODUCTION 

o 

Ferrofluids ||] are oil- or water-based colloidal suspensions of permanently magnetized particles. In an applied 
J3 magnetic field the particles align creating a strong paramagnetic response in the ferrofiuid. Because they are flu- 
ids, these siispensions can flow in response to forces. For example, ferrofluid droplets elongate parallel to applied 
fields and undergo tip-sharpening transitions . When a ferrofiuid droplet is confined between two plates 

^ \ in a "thin film" geometry, surrounded by an immiscible fluid, and a fleld is applied perpendicular to the plates, 
Q it undergoes fleld induced bifurcations [|l0| leading to intricate labyrinthine patterns Ferrofluid emulsions 

undergo structural transitions under an applied field from a randomly dispersed structure of the emulsion droplets 
to droplet chains, columns and worm like structures depending on volume fraction, sample geometry and the 

rate of field application. 

^ , A droplet of ferrofluid elongates under applied field because of the demagnetizing flelds of magnetic poles on the 
^\ • surface of the droplet. Surface poles arise wherever the droplet magnetization has a component perpendicular to the 
OO \ surface. The demagnetizing fleld that they create opposes the magnetization, creating a demagnetizing energy that 
^""^ ■ depends on the shape of the droplet. The droplet elongates to reduce its demagnetizing fleld and energy. Because 

' elongation increases the surface energy of the system, an equilibrium shape is reached when the magnetic forces 

, balance against the surface tension forces. 

' The elongation of freely suspended, 3-dimensional droplets has been well studied [||-|5| . The droplets can be assumed 
, to be ellipsoids for small elongation. The demagnetizing fleld is thus uniform and the elongation (major axis minus 
minor axis divided by minor axis) is found to be proportional to the undeformed droplet radius. The case of droplets 
^ conflned in "thin fllm geometry" however, involves two length scales, droplet thickness and its undeformed diameter. 
I In the limit of small aspect ratio (droplet thickness divided by its undeformed diameter) the demagnetizing fields are 
■ O stronger near the edges of the droplet than at its center. We flnd that the elongation divided by droplet thickness in 
^ . this geometry is proportional to the logarithm of the aspect ratio. 

Prior experiments ||^,^ have proposed droplet elongation as a tool for measuring surface tension between the 
ferrofluid and the surrounding immiscible fluid. We improve on the existing theory ^ by incorporating spatial 
variation of the demagnetizing fleld inside the droplet. We perform an experiment supporting our predicted logarithmic 
behavior. 

Section |l^ of this paper presents our theoretical study of the elongation of a ferrofluid droplet conflned within a thin 



o 
o 



X 



■ film. Our principal result is a predicted logarithmic dependence of elongation on droplet aspect ratio. We contrast 



this result with the corresponding elongation of unconfined droplets. Section III describes an experiment done with 



ferrofluid emulsions that tests our theory. The experiment is in qualitative agreement with our theoretical prediction. 



but differs quantitatively in at least one respect. In section IV we discuss a possible explanation of the discrepancy 
based upon droplet contact angles with the conflning plates. 
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II. THEORY 



Consider a paramagnetic liquid droplet confined in a thin film between two parallel plates with a gap, A, in the z 
direction (see figure 1). An immiscible liquid surrounds the droplet. Let the thickness, A, be much smaller than the 
radius of the undeformed droplet, tq. This small aspect ratio 

f-f (1) 

provides the pseudo-two-dimensional character of the problem. If a uniform, weak, field Hq is applied parallel to the 
plate, the droplet magnetizes. The magnetization creates an opposing demagnetizing field whose strength depends 
on the droplet shape. The droplet elongates to decrease its magnetic energy, reaching equilibrium when the magnetic 
forces balance against the restoring forces due to surface tension. In this section we define the elongation of the droplet 
and calculate the surface energy, i?5, and the magnetic energy, E'a/ , of the droplet as a function of its elongation. By 
minimizing the total energy with respect to the elongation we obtain the elongation as a function of Hq, ro, and A. 

For simplicity assume the elongated droplet has a uniform cross section, C, independent of z. This corresponds to a 
contact angle of 90° between the paramagnetic liquid, the surrounding fluid and the glass plates, and a plate spacing 
much less than the capillary length of the two liquids. Thus the droplet has straight edges if viewed from the side 
(see figure 1). The role of contact angle will be discussed later in section |^ We write the equation for C in polar 
coordinates as a generic smooth perturbation to a circle, 

r = ai + Q!2 cos 26*. (2) 

We only include a single harmonic, since we expect coefficients for the higher harmonics to be much smaller than 
a2 for small perturbations. The cross section C has semi- major axis a, and semi-minor axis h (see figure lb), with 
ai = (a -I- 6)/2 and a2 — {a — b)/2. We define the elongation of the droplet 

e.--l. (3) 

We assume that the elongation, e, is much less than 1. Imposing the constraint that the volume of the droplet (A 
times cross-sectional area) remains constant we calculate 

(l + F/2)l/2' "2 (1 + ^2/2)1/2' W 

where k = e/{2 + e). 

The surface energy is the sum of interfacial areas times surface tensions between all pairs of the three phases (solid 
glass, ferrofluid droplet and immiscible fluid). For the case of uniform cross-section (90° contact angle) droplets, the 
glass-ferrofluid and glass-immiscible fluid interfacial areas are independent of the shape of C due to the fixed volume 
constraint. Hence we concern ourselves with the droplet-surfactant solution interface, the area of which is A times 
the perimeter. The perimeter of cross section C can be calculated as a power series in e, 

5 = 2^ro(l + ^62 + 0(e3)). (5) 

As expected, the leading correction to S is second order in e since the perimeter should increase regardless of the sign 
of e. The relevant surface energy of the droplet is 

Es = apiSA (6) 

where api is the surface tension of the ferrofluid-immiscible fluid interface. 

The total magnetic energy of any paramagnetic body under applied field is [ p^ 

Em = -lJ^d^rUo -Mir). (7) 

The magnetization M(r) is determined by the self consistent equation 

M(r) =x(Ho+Hz,(r)) (8) 
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for linear susceptibility x, where 

H^(r) = 1^ £v' (M(r') • n{v'))^^^ + (V • M(r'))^^-^ (9) 

is the demagnetizing field due to the magnetization M(r), with n(r') being the outward normal at any point on 
the surface. The surface integral gives the demagnetizing field due to the surface poles which appear wherever the 
magnetization has a component normal to the surface. The volume integral gives the contribution to the demagnetizing 
field due to volume charges which appear at points where the magnetization has non-zero divergence. 
To calculate the magnetic energy we expand M and Hd in power series in the susceptibility X; 

M(r) = M(i)(r) + M(2)(r) + M.^^\v) + ... (10) 



Un{v) = HW(r) + Hg)(r) + Hg)(r) + (11) 



where M'^")(r) and (r) are proportional to x"- Equating terms in (g) of equal order in x we get 

M(i)(r)=xHo (12) 

and 

M("+i)(r) = xH(3"\r). (13) 

Note that M*^^^(r) is independent of r because the applied field is uniform whereas M^"''(r) may depend on (r) for 
71 > 1 because Hd (r) may be non-uniform. To second order in x we write the magnetic energy of the droplet in (Q) 
as 

Em = -IxH'oV~1 |^d^rMW.H«. (14) 

The first term in equation ( p^ ) for the magnetic energy is independent of the shape of the droplet and hence 
unimportant for our consideration. The second term in the energy is the demagnetizing energy Ejj due to a uniform 
magnetization M^^^ = xHq- Because M*^^) is uniform there are no volume charges, and the surface poles appear only 
along the droplet-immiscible fluid interface, to first order in x- Rewrite the second term in ( p^ as an energy due to 
the induced surface charges along the curved surface of the droplet 



2 , /""^ ,//,/,/ (n • Ho)(n' ■ Ho) 
" 2^ J J J J |rr^| 



(15) 



Here ds and ds' are infinitesimal arc- lengths along the contour of the droplet C, and n and n' are the outward normals 
to the curved surface of the droplet at points {s,z) and (s',z') respectively. 

Write |r — r'l = R"^ + {z — z')^, where R is the in-plane distance between points at positions s and s' on C. 
Integrating over z and z' in ( [Tsl ) yields ||l^ 

Ed = x^A <fds(j^ ds'{h ■ Ho)(n' • Ho)$(i?/A) (16) 

where 



$(i?/A) = R/A - + (i?/A)2 + In [{R/A)/{^/iT{R/W - 1)J • (17) 

Using equation (^) for C we calculate the demagnetizing energy in ( p^ ) as a series expansion in e and the aspect ratio 
P = A/2ro 



ED=X^H^v\2p\n--3ep\n- + ---] (18) 
I p p ) 



where V = Trrg A is the volume of the droplet, and B = Ae and C — Ae are geometrical constants. The term 
in the brackets can be identified as 27r times the demagnetizing factor of the droplet along the direction of applied 
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field. Additional terms in the series in equation ( |i8| ) are of higher order in e or in p. For small elongation and large 
aspect ratio we may neglect these higher order terms. 

Minimizing the total energy E — Es + Em with respect to e gives 

e=^ln^. (19) 
api p 

Corrections to this result are higher order in aspect ratio p or higher order in e itself. Interestingly, the elongation 
depends only logarithmically on the undeformed radius tq, and has a much stronger dependence on the thickness, A, 
of the droplet. This result differs from an earlier theory |g] which omits the logarithm because it assumes that the 
demagnetizing field is uniform inside the droplet. 

In the case of unconfined, nearly ellipsoidal droplets the demagnetizing field is quite uniform inside the droplet. 
The demagnetizing energy is therefore proportional to the volume ((4/3)7rrQ) of the droplet according to equation (pH). 
The surface energy is proportional to the area (47rrQ) and the elongation is thus proportional to tq. In the case of thin 
film geometry, however, the demagnetizing field is very non-uniform. For distances much less than A near the droplet 
edge, the component of the demagnetizing field is of order M, since the edge acts like an infinite sheet of charge in the 
first approximation. For distances much greater than A the demagnetizing field is of order MA / r since the edge acts 
as a line charge in this case. The contribution to the integral for the demagnetizing energy in equation (|l^) mainly 
comes from the bulk of the droplet and goes like TqA^ ln(ro/A). The surface energy is proportional to 27rroA and 
the elongation is therefore proportional to Aln(ro/A). The logarithmic variation of elongation with the aspect ratio 
is thus a signature of the non- uniform nature of the demagnetizing field inside the droplet. 



III. EXPERIMENT 



A. Setup 

a. Sample Preparation and Structure Our sample consisted of a ferrofluid/ aqueous solution emulsion confined 
between two glass plates. The oil-based ferrofluid used was EMG 905 made by Ferrofluidics. To reduce the surface 
tension between the ferrofluid and the immiscible aqueous external phase, we incorporated surfactants in the aqueous 
phase. A solution of a commercial detergent made the best emulsions while solutions with other pure anionic surfac- 
tants either showed hardly any elongation of the ferrofluid droplets under applied fleld or produced droplets without 
sharp boundaries with the aqueous phase. In contrast, our stable, well behaved emulsions allowed us to probe and 
confirm the fundamental aspects of our model. 

To prepare the emulsions, a single drop of ferrofiuid (~ 0.1 ml) was added to 10 ml of surfactant solution which was 
a 12 times dilution of the commercial detergent. The liquid was shaken (by hand) to prepare the emulsion, creating 
ferrofiuid droplets with diameters varying from ~ 5 — 200 /im. A small amount of this emulsion was then put between 
two glass plates which were circular, about 2 cm in diameter and 4 mm in thickness. These plates were cleaned using 
soap and alcohol and then rinsed with ROPure water. We also tried acid cleaning of the glass plates, however it did 
not result in any noticeable change in the quality of the sample. 

We used a rectangular spacer made of mylar foil to separate the plates and prevent the emulsion from leaking 
out from the edges of the plates. The mylar foil extended to to the edges of the glass plates and had a rectangular 
hole in the center into which the emulsion was inserted. The thickness of a single mylar spacer was measured to be 
6.54 ± 0.06 /im. The experiment was performed with one and two spacers to ensure small aspect ratio. 

For the cell assembly, the mylar spacers were placed on the first plate and a drop of the emulsion was put in the 
center of the plate. The second plate was placed on top and the two plates were clamped together using a pair of 
brass rings. The rings were tightened by a set of 4 equally spaced screws. We measured the thickness variation across 
the sample by making a "dry" sample (without the emulsion) and counting resulting white light interference fringes. 
Although the thickness of mylar spacers was measured to an accuracy of 1 percent, the thickness variation across the 
sample was found to be 10% resulting from the stresses due to clamping and possible entrapment of dust in the cell. 

b. Apparatus A schematic diagram of the experimental setup is shown in figure 2. We put the sample at the center 
of a pair of Helmholtz coils to insure a homogeneous magnetic field. The field measured close to the sample using 
a Hall probe showed a variation of less than 4% across the sample. The sample was set up horizontally to prevent 
gravitational settling of the ferrofiuid droplets. Horizontal alignment was achieved using a bubble level. 

The sample was illuminated from below using a diffused light source and observed from above using a tele- 
microscope. The tele-microscope was connected to a CCD camera and the image from it was fed into a video 
recorder and recorded on video tape. Images from the recording were later processed using NIH Image. We calibrated 
the optical system using a measuring reticule aligned along the two orthogonal directions of the CCD array. Figure 
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3 shows a low magnification view of a typical sample. The ferrofluid droplets appear much darker in the image than 
the surfactant solution around them. 

c. Experimental Procedure and Image Analysis During the experiment the applied field was incremented every few 
seconds. We found the response of the droplets to the field to be nearly instantaneous and the shape of the droplets 
remained constant at constant field. Experiments with decreasing field strength showed no hysteresis in droplet shape. 
While droplet elongations were observed to be small we incremented the field in steps of about 1 Gauss, and increased 
the increments up to about 5 Gauss as the elongation increased. Droplet elongations appeared to vary smoothly with 
applied fields over the entire range from to 50 Gauss. 

During each experiment the droplets were observed on a video monitor and recorded on tape. After grabbing images 
of distorted droplets, we used a cut-off in pixel gray scale level to identify the droplet edge. The semi-major axis (a) 
and the semi-minor axis (6) were directly read off the image using NIH Image. At zero field measured elongations were 
small (RMS magnitude around 0.003) and in random directions. These minor perturbations from a circular shape 
were likely due to microscopic distortion of the contact line pinned on weak surface heterogeneities. The "observed 
radius" rg was calculated as the average of the two semi-axes at zero field and the elongation at each field value was 
calculated using data analysis software. 

d. Results For each of the 48 droplets studied we plotted elongation e, versus the square of the applied field, Hq. 
Figure 4 shows typical plots. The elongation is proportional to the square of the applied field for small applied fields 
as predicted. Saturation effects, although small, can be seen at higher values of the field. The plot of elongation, for 
each droplet was fitted to 

^ = ko + kiHl + k^Hl (20) 

We included terms only up to order because the saturation effects were observed to be small. We include to 
allow for the observed small elongations at zero field. 

The coefficients ki of each droplet were then plotted versus the inverse of the aspect ratio 1/p ~ 2ro/A on a 
semi-log plot (see figure 5). The theory predicts a slope of /cffi and an intercept of 1/C on the horizontal axis 
with C = 1.74. The data points in figure 5(a) fall on a straight line as predicted by the theory. Also, as predicted 
by the theory, the data points for two different droplet thicknesses overlay each other. There is substantial scatter 
in the data, but the deviations from a straight line are random and consistent with the error bars. The chief source 
of uncertainty was the 10% uncertainty in thickness due to the variation observed across the sample. Figure 5(b) 
displays the deviation of ki from the best fit normalized by the uncertainty. The uncertainties in measuring e, ro and 
Ho were found to be negligible in comparison. 

Dividing the susceptibility x = 1-9 for the ferrofluid used |jl^ by the slope = 0.119 ± 0.004 cm/dyne obtained from 
the fitted line we get api = 30.4 ± 1.1 dynes/cm, typical of oil-water surface tensions. From the fitted line we also get 
C = 0.35 ± 0.08, differing substantially from our theoretically predicted value of 1.74. It may be possible to explain 
this discrepancy by considering the effect of the contact angle of the ferrofluid-immiscible fluid interface with the glass 
plates. In the discussion section below we explore the qualitative effect of the contact angle. 

In figure 6(a) we plot e/A versus 2ro/A on a linear scale. If the demagnetizing field inside the droplet was uniform 
like in the case of unconfined droplets, the plot would be a straight line. However, the plot is clearly not a straight 
line and the deviations from the best fitted straight line are systematic (see figure 6(b)). This further supports our 
theoretical result that the demagnetizing field inside a confined droplet is non-uniform and the elongation divided by 
thickness is proportional to the logarithm of the aspect ratio. 



IV. DISCUSSION 

The results discussed in section [II agree with our theoretical prediction ( p^ ) of logarithmic variation of e/A with a 
proportionality constant of x^ japi- However, our theoretical value for C is 4e^^/^ = 1.74 whereas the experimentally 
measured value for C is 0.35 ± 0.08. 

One possible explanation for the discrepancy in the value of C is that the ferrofluid/glass contact angle is not 90° 
and consequently the cross-section of the droplet is not uniform. Our calculations are for uniform droplet cross-section, 
which corresponds to a contact angle /? = 90° between the glass plate and liquid droplet. The experiment, however, 
was performed with an oil-based ferrofluid in a surfactant solution for which the oil-glass contact angle < 90° (see 
figure 7). A contact angle of other than 90° will affect the elongation in two ways: by changing interfacial areas to 
alter the functional form of Es and by redistributing the magnetic surface poles to alter the functional form of Em- 
We consider these two effects in turn. First, however, we must address an ambiguity in the definition of aspect ratio 
and elongation which results from the non-uniformity of droplet cross-section. 
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Our experiment observes the profile of the largest cross-section of the droplet. For a circular droplet with (3 < 90° 
this is the radius ri defined as the radius at mid-gap as shown in figure 7. For an elongated droplet we measure the 
semi- major and -minor axes ai and bi and, through equation (||), the elongation ei. We also define r2, 02, 62 and 62 
associated with the ferrofluid- immiscible fluid-glass plate contact line (see figure 7). Since A is much less than the 
capillary length of the ferrofluid/immisciblc fluid, to a good approximation |0 the profile of the droplet will be an 
arc of a circle, so the difference between ri and r2 is of order A, and likewise for the semi-major and -minor axes. 
The difference ei — £2 is of order A/ri relative to the elongation. Recall that our result ( p^ ) for the elongation is only 
the lowest order term in a series expansion in the aspect ratio. Thus the distinction between ri and r2, and between 
ei and £2, does not alter our result at the lowest order in aspect ratio. 

When the contact angle differs from 90°, the cross section of the droplet depends on z. Consequently, the contact 
areas of the glass plates with the droplet and with the surfactant solution may vary as the droplet elongates. All the 
three interfacial areas must be taken into account to calculate the surface energy. The total surface energy is 

Es ^aFiAc + 2(jFGAG + '2(JiG{A- Ag), (21) 

where the three surface tensions between ferrofluid-immiscible fluid, ferrofluid-glass, and surfactant solution-glass, are 
denoted by api, <yFG, a-nd ct/g respectively, Ac and Ag are defined below and the total area of the sample is denoted 
by A. The factors of 2 in the second and third terms of the surface energy account for the two glass surfaces. 

The area of the droplet-surfactant solution interface Ac is given approximately by the circumference of C multiplied 
by the arc length of the bulge 

Ac = 2..,(l + Ae^)Ai!^. (22) 
16 cosp 

We use ri here to calculate the circumference of the droplet because is the radius observed during the experiment. 
To first order in the aspect ratio, using ri or r2 in equation (|2|) yields the same result. 

The droplet's contact area with the glass plates must be adjusted to maintain a constant total volume of ferrofluid 
as the droplet elongates. We approximate the volume of the bulging region by the circumference of C multiplied by 
the projected area of the bulge. The contact area Ag must be adjusted so that AqA changes by the negative of the 
change in volume of the bulge. Thus we write 



Ag = 27rr2 



3 f (7r/2-/3) ^\ A 2' 



1 ~ ^\ Y7i tan/3 ^— e 

32 cos"^ p J r2 



(23) 



Using r2 instead of ri makes the above result exact for zero elongation. The area of ferrofluid in contact with the 
glass plates decreases with elongation for an acute contact angle because the volume of the fluid contained in the 
outward bulge of the droplet increases and therefore the fluid contained in the bulk of the droplet decreases. For 
obtuse contact angles exactly the opposite happens for similar reasons. 

To understand how the contact angle affects the magnetic energy, consider the work done by the magnetic field 
as we change the contact angle from 90° to (5 while keeping the volume of the droplet constant. This work, divided 
by the circumference , must be independent of ri in the limit ri going to infinity, since the magnetic field near the 
surface of the droplet will not depend on ri in the large ri limit. The work done by the magnetic field is the difference 
in energy between the straight edge droplet with contact angle of 90° and the bulging droplet with a contact angle 
of /3. The demagnetizing energy of the bulging droplet must therefore have the same dependence on ln(ri/A) as the 
straight edge droplet or the difference in the demagnetizing energies divided by the circumference will be proportional 
to ln(ri/A) and will blow up in the large ri limit. Hence, the demagnetizing energy for the bulging droplet must be 
identical to equation (|l8|) but with different values B and C replacing the constants B and C. 

As the droplet bulges inward or outward the charges on the surface get distributed over a larger area, decreasing 
the demagnetizing energy. The constant B therefore has a smaller value for a bulging (inward or outward) droplet 
than _B, the value for a straight-edged droplet. However, since the demagnetizing energy is always positive, smaller 
demagnetizing energy {B < B) implies a weaker dependence of demagnetizing energy on elongation. Thus, we expect 
the value of C to be smaller for a bulging droplet than the value C for a straight-edged droplet. 

Finally, consider how the contact angle dependence of surface and magnetic energies affect the elongation calculated 
in equation (|l^) for the case (3 = 90°. The e dependence of the surface energy remains quadratic, but the coefficient 
now depends upon a linear combination of the three surface tensions ctf/, (Jfg ^md ct/g- This combination will replace 
a PI in equation (p^. The functional form of the magnetic energy remains unchanged, but the values of B and C 
depend on contact angle. Thus the smaller value C replaces C in equation (p^). For (3 7^ 90° the experiment cannot 
be used to determine aps unless ctfg and uig are known. Since the (3 in general is not 90°, it is only possible to 
measure the effective surface tension during elongation, and not api itself. 
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V. CONCLUSIONS 



We study the elongation ferrofluid droplets, confined in thin film geometry, under weak applied field. Our theoretical 

calculations predict the elongation of a droplet depends logarithmically on aspect ratio. This behavior contrasts 
with the case of unconfined 3-dimensional droplets where elongation is directly proportional to undeformed droplet 
radius. We measured the elongation of ferrofluid droplets in an experiment performed on ferrofluid droplets in a 
ferrofluid/water/surfactant emulsion. The results of our experiment agree with the functional form our theoretical 
prediction, however the experimentally measured value of C differs from the predicted value. We suggest the droplet 
contact angle with the confining plates as a source of this discrepancy. 
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FIG. 1. (a) A side view of a ferrofluid droplet confined between two glass plates, (b) A top view of a ferrofluid droplet 
elongating under applied field. The dashed line shows the undeformed droplet. N and S indicate the north and south magnetic 
poles. Ho is the demagnetizing field. 



CCD 



TM 



VR 



HC 



HC 



LS 



FIG. 2. A schematic diagram of the experimental setup. LS=Light Source, HC=Helmholtz Coil, S=Sample, 
TM=Tele-Microscope, CCD=CCD Camera, VR= Video Recorder, C=Computer 
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FIG. 3. A small magnification view of the sample showing ferrofluid droplets in emulsion. 
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FIG. 4. The plot of elongation vs for droplets with different radii and two different thicknesses. The error bars are 
smaller than the size of the symbols on the plot. The symbols and radii for droplets with thickness spacing of A = 6.5 /im are 
(x, 96.5 nm), (•, 112.0 /xm) and (+, 216.5 /xm). The symbols and radii for droplets with thickness spacing of A = 13.1 iim 
are (□, 50.0 /um), (■, 98.0 fim) and (©, 177.0 /um). 
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FIG. 5. (a) The plot of fei vs 1/p = (2ro/A) on a log scale. The dashed line is the best fit of the data to a straight line, (b) 
The deviation plot of the data from the best fit normalized by uncertainty of each data point. 
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FIG. 6. (a) The plot of fei vs 1/p = 2ro/A and the best fit straight Une. (b) The deviation plot of the data from the best 
linear fit. 
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FIG. 7. A ferrofluid droplet making an acute contact angle with the glass plates. 
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